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At finite temperature and in non-equilibrium environments we have to resum perturbation theory 
to avoid infrared divergences. Since resummation shuffles the perturbative orders, renormalizability 
is a nontrivial issue. In this paper we demonstrate that one can modify any type of resummations 
- even if it is momentum dependent, or it resums higher point functions - in a way that it is 
renormalizable with the usual zero temperature counterterms. To achieve this goal we reformulate 
resummation in form of a renormalization scheme. We apply this technique to perform a 2PI 
resummation in the six dimensional cubic scalar model in equilibrium. 
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I. INTRODUCTION 
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Renormalized perturbation theory in MS scheme was developed to extract perturbative information about vacuum 
expectation values of observables at high energies. Still, when it fails, for example in computation of expectation values 
with a generic density matrix, we tend to speak about "non-perturbativc effects" and the failure of perturbation theory 
£NJ ■ (PT). An example at high temperature when this pathological behavior occurs was demonstrated by Dolan and Jackiw 
^ ! Hi showing that in <E> 4 theory at finite temperatures higher order diagrams can provide contributions ~ (AT 2 /m 2 )™ 
which is large at high temperatures. They also found the solution of this problem: with (re)summation of a specific set 
of diagrams (the daisy diagrams) , the problematic ~ XT 2 term becomes part of the propagator mass, thus appearing 
. in the denominator instead of the numerator. 

This solution reveals, at the same time, the core of the problem. Perturbation theory in MS scheme starts with 
\Q , a free action, fitted to the zero temperature particle spectrum. At finite temperatures the true excitation spectrum 
changes significantly. Since particle spectrum comes from the free part of the Lagrangian, any perturbation theory will 
fail which tries to follow this change via interactions. This line of thought also provides an easier way of performing 
daisy resummation, by changing particle spectrum of the theory from the beginning. Technically one adds to and 
subtracts from the Lagrangian the same (temperature dependent) mass term, treating one as a part of the free action, 
the other as a (thermal) counterterm [2j . The modification mass term can be set equal to the high temperature limit 
of the self energy diagram, or it can be optimized via a gap equation 0-01 

Taking into account thermal effects only through an effective mass can be used only cautiously in case of more 
complicated situations, as in non-equilibrium, or in gauge theories |21]. In these cases the non-trivial, non-quadratic 
momentum dependence of the self energy is important, too, and it must be taken into account already in the free 
Lagrangian. In gauge theories the hard thermal loop (HTL) improved perturbation theory [f| provides this type of 
resummation, which uses the high temperature limit of the momentum dependent self energy as part of the effective 
propagator (and uses also effective vertices to maintain gauge invariance, cf. [j|). In the generic case the modification 
of the spectrum can be taken into account most effectively with the method of Cornwall, Jackiw and Tomboulis Q, 
which performs two-particle-irreducible (2PI) resummation. In this method we should compute the diagrams with 
the exact (also background field dependent) propagator, and, consequently, no self-energy corrections should be taken 
into account. This method was intensively developed and used in numerical calculations in the recent literature 0. 

The consequent application of all these methods is obstructed by their non-trivial relation to the renormalization. 
Direct application of the thermal counterterm method, for example, requires temperature dependent infinite countert- 
erms to cancel the one-loop level ~ m^/e type of divergences, where m\ is the thermal mass 0, 8, 9]. The temperature 
dependent piece drops out at two loop level, but there new temperature dependent divergences are generated. The 
situation is even more serious in the case of 2PI method: at one loop level we generate terms to the exact propagator 
that are ~ m 2 logp; when putting back into loops, this term provides us with new types of divergences, not present 
in the original renormalized Lagrangian. Therefore the resummation obtained by truncating the 2PI equations at a 
given loop order is not consistent as a renormalizable theory. 
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The main reason of failure is that we separate the loop order of the IR sensitive diagram we want to resum and 
the corresponding counterterm diagrams, which, in the original PT, made them finite. In a resummation based on 
fixed order skeleton diagrams this leads to uncompensated divergences. There are different methods to keep track 
of the divergences. The local case is simpler and there we have a deeper understanding of the problem. A self- 
consistent method for resummation was developed in |8l Il0| to obtain finite gap equations. In 11] the self-consistent 
resummation in the framework of 2PPI effective action is accomplished by explicitely introduced counterterms into 
the Lagrangian. In H this method is generalized to any, momentum-independent resummations. Here an alternative 
viewpoint is applied: the resummation being treated as a change of renormalization scheme. The idea of adapting 
the renormalization process to the environment was first suggested in |12| . 

The momentum dependent resummations are more difficult, here the results are concentrated on the 2PI resumma- 
tion. In a seminal paper, van Hees and Knoll ^| presented a proof based on the BPHZ scheme showing that it is in 
fact enough to use the zero temperature BPHZ scheme to renormalize the finite temperature 2PI resummation. The 
central problem there is to renormalize the zero tem per ature ladder resummation, ie. the Bethe-Salpeter equations. 
This method was further developed and used in [Tj, [l]| [T^| . A method to renormalize Hartree-Fock approximation 
was developed in |l7j . 




It is, however, still an open question, what kind of momentum dependent resummations are renormalizable. In this 
paper we study this question, primarily in the framework of the four dimensional $ 4 and six dimensional < & 3 model, 
but the logic can be easily extended to more complicated models. The starting point is the method of 9], which is 
generalized to momentum dependent case. The main idea is that, instead of resumming the necessary counterterm 
diagrams in the original resummation, we slightly change the resummation procedure: in the low energy (IR) regime 
it will be still the desired resummation, while at high energies (UV) it falls back to the ordinary perturbation theory. 
This means that the divergences remain identical to the usual zero temperature ones. 

The paper is structured as follows. In Section ^ we give a different view of resummation, interpreting it as 
a specific renormalization scheme even in the momentum dependent case. We generalize somewhat the notion of 
"renormalized Lagrangian" , allowing momentum dependent parameters, and taking care not to spoil the perturbative 
renormalizability. In Section IlIII we summarize the results of 2PI resummation, and find that specific resummation 
scheme that corresponds to it. It will be shown that it is not a well-defined scheme (as it can be conjectured from the 
renormalization problems of the 2PI resummation itself), so in Section lill Bl we suggest improvements of the method 
that is capable to satisfy renormalizability and 2PI resummation in the IR regime at the same time. In Section II VI 
we apply the method to find the 2PI resummed self energy of the six dimensional <I >3 model. In Section we give 
our concluding remarks and discuss open questions. 



Resummation is a notion that pre- assumes a singled out perturbation theory (the "unresummed" one). It suggests 
that in this singled out perturbation theory (PT) we have to do additional work to compute some observables by 
taking into account contributions from higher orders. Usually this singled out PT is the one stemming from the MS 
renormalized Lagrangian. This works nicely when used for calculation of vacuum expectation values of high energy 
observables, but needs to be resummed, for example, at finite temperatures. 

This way of thinking is, however, misleading, as we have no a priori reason to prefer one environment/process 
above others. Having in mind this democracy of environments, we should not speak about resummation, but about 
perturbation theories best adapted to a certain environment (cf. also |T^L 

The definition of quantum field theory is based on the path integral living on a fine enough spacetime lattice, 
with an action (the bare action) that have parameters depending on the lattice spacing. This UV action provides, in 
principle, the expectation values of any IR observables, at any environment. If we do PT, however, we must have some 
knowledge about the result in the specific environment in order that already the first approximation is "sufficiently" 
close to the exact result: ie. we should ensure that the corrections are indeed "perturbative" , smaller than the leading 
order. Since in the usual PT the solvable part is the free (quadratic) theory, the above requirement says that the 
vacuum as well as the free excitation spectrum should be sufficiently close to the exact one. But the vacuum of the 
theory may change in different phases of the matter, and also the excitation spectrum can be dramatically different 
at different physical environments, for example at different temperatures, not to speak about the non-equilibrium 
environment. So, although it is implicit in the usual discussions, before we start PT, half of the task should be solved: 
to guess the vacuum and the spectrum sufficiently well. 

Different PTs, adapted to different environments, yield different results already at tree level. But, of course, at 
infinite order any consistent PT should provide the same result for a given observable. That means that the result 
of the tree level of PT A is somehow hidden in the higher orders of PT B. From the point of view of B, therefore, 
A represents a resummation of the perturbative series - and vice versa. And both A and B are resummations of the 
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bare perturbative series. A zero temperature example is that results obtained in the on-mass shell scheme seem to be 
"resummed" from the point of view of the MS scheme. 

In addition to the necessity of resummation there is another problem. When we start with a UV theory with 
appropriate bare parameters, it is granted that any IR observables are "finite", ie. do not depend on the exact value 
of the UV cutoff. The requirement for a consistent PT is, however, much more severe: we should have finite results 
order by order. Only a very smart resummation of the bare theory known as the renormalizcd PT, can provide finite 
results at each perturbative order |Tsj . No other perturbative methods are known that would be capable to render 
the theory IR finite. Therefore, if renormalizability is an issue, we must not use other PTs than the renormalized 
ones. 

Technically, in renormalized PT, we split the parameters into a renormalized and counterterm part, the latter being 
cutoff dependent, and should formally be taken into account at higher orders of the calculation. Renormalizability 
requires specific form only for the UV part. Thus it allows to work with different PTs, which use different finite parts 
of the counterterms: these are called renormalization schemes (RSs). When we want to fit the PT to the external 
environment the only thing we can play with is the finite part of the counterterms; ie. we should adapt the RS to the 
environment by choosing appropriate finite parts. 



A. Class of constant physics 

In the family of RSs we can define an equivalence relation: two RSs are equivalent if they provide the same 
expectation values for any n-point functions at infinite order (up to a finite normalization factor). This leads to the 
equivalence classes of schemes that provides "constant physics" . 

Let us take two RSs: A with parameters gf, and B with parameters gf . The expectation value of an n-point 
function O n = T$(xi) . . . &(x n ) will depend on these parameters: 

(d n ) A = A (g A ), {6 n ) B = B {g B ). (1) 

A and B belong to the same class of constant physics, if we have 

A ( 9 A )=a B B (g B ), (2) 

with a (n-independent) constant Cab- According to the theory of renormalization |Tsj these equations can be fulfilled 
for any n, if there is a specific relation between the renormalized parameters 

gf = f t AB (9 B )- (3) 

Sometimes these equations are referred to as the renormalization group (RG) equations. The physical condition for 
providing the same values for the observables from both schemes is that they both stem from the same bare Lagrangian 

To demonstrate these notions we give an example in the $ model. The Lagrangian reads in terms of the bare 
quantities (denoted by subscript "0"): 

C = \{d^?~^l-^i (4) 
The renormalized Lagrangians are used to write in the following form 

£=-(<9$) $ $ H (<9$) $ $ , (5) 

2 V ' 2 24 2 V ' 2 24 ' K ' 

where the formal power of A in the counterterms determine the loop order they first contribute. Since determines, 
nonperturbatively, the expectation values of all IR observables, clearly the renormalized Lagrangians of J5J are giving 
the same value (at infinite order) , if we satisfy 

$ = Z$, Z 2 ml = m 2 +6m 2 , Z 4 X = X + 5X, (6) 

where Z 2 = 1 + SZ 2 . Therefore two schemes give the same result (up to normalization), if they obey: 

Z a = CabZ b , m A + 6m A = ( AB (m 2 B + 6m 2 B ), \ A + 5\ A = ( AB (\ B + 6\ B )- (7) 
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The condition that schemes A and B belong to the same class of constant physics can be rewritten compactly, using 
the generator of 1PI diagrams as 

T A [§\=T B [U B §]. (8) 

At finite nth order this is satisfied up to order g n . Then the left hand side contains contributions from higher 
orders, as considered from RS B. So the left hand side provides a resummation for RS B. Moreover, since all the 
functions f AB and A are finite, we have a finite, ie. renormalized resummation. 



B. Momentum independent resummations 

Let us try to apply the ideas of the previous subsection to the momentum independent mass resummation in <J> 4 
model. A more comprehensive study of this case can be found in Ref. Q . 

Let us choose for RS A the MS scheme. At one loop level with dimensional regularization the counterterms are 
defined as M : 



SZ MS,l = °> 



A 



32tt 2 



— h 7b — 1 — ln47r 



MS,1 



3A 2 / 1 , , 



In this scheme zero temperature observables can be calculated with reasonable precision; at finite temperatures, 
however, IR divergences make the calculation difficult. For example the 1-loop self energy at zero momentum reads, 
in dimensional regularization: 



n MS (P = °) = ™ 2 



Am 
32^ 



2 r 



1 m 2 
— + 7e - 1 + m l — o 



A 



J(m) + Smfjg , 



where /i is the renormalization scale. J is the finite temperature part of the bosonic tadpole function: 



(10) 



J(m) 



2tt 2 



dui V 1 u 2 — to 2 n(oj), 



(11) 



where n(uj) is the Bose-Einstein distribution. It can be approximated at high temperatures as 



J(m) * 12 ~ 



1 9 

m m 

Te^ 2 ^df 2 



where lnc = 1 + 21n47r — 27^ ps 4.90762. Therefore, at high temperatures, the one-loop self energy will be 

XTm 



n M s( P = o,T) 



XT 2 
~24~ 



Am 2 , cT 2 
■ In ■ 



(12) 



(13) 



MSV1J "" ' 24 8tt ' 32tt 2 *" Airfi 2 ' 

This means that the first order correction can be larger than the tree level term, if AT 2 3> m 2 , ie. at high temperatures. 
At higher orders this failure of the MS scheme becomes more and more obvious, as we receive (AT 2 /m 2 )™ contributions 
from the daisy diagrams 0. 

This phenomenon, however, does not mean the failure of perturbation theory itself, since we can easily find a RS, 
where this problem does not show up. Indeed, the finite temperature version of the on-mass-shell (M) scheme uses 
the counterterm 



Sm 2 M1 



Am 
"32^ 



2 r 



IE 



1 + In ■ 



47T^ 2 



^J(m), 



and so the self energy in this scheme reads simply 

n M (p = o,T) 



(14) 



(15) 



This is just the tree level value - as we expect from the on-mass-shell scheme. No dangerous AT 2 terms appear, and 
by definition it extends to all orders in PT. So it seems to be much better to use the on-mass-shell scheme at high 
temperatures for perturbative calculations. 
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If we want to compare the results of the MS and the finite temperature on-mass-shell schemes - for example if we 
want to follow the temperature dependence of the quantities starting from zero temperature - then we must take care 
that they belong to the same class of constant physics. Since now only the mass counterterm is changed, we have 
C = 1 and Aj^g = Am = A, and we must satisfy 



_ Am Ms 
M s 32tt 2 



1 



7b — 1 — ln47r 



C(A 2 



32tt 2 



1 , TO' 

— +7B-l + ln-^ 
e An^ 



A 



At 0(A) this simply yields 



32tt 2 



J(tom)- 



J(to m ) + 0(A 2 ). (16) 



(17) 



When we solve this equation - a finite gap equation - the result will be applicable in the given class of constant 
physics. In particular, we can interpret this result as a resummed MS self energy: 



yrresummed 
MS 



(P = 0,T) = rn 2 M (to^jo, A) 



(18) 



An alternative to this approach is when we define a scheme where only the dangerous finite temperature part is 
resummed. This F scheme is defined by the counterterm 



5m 2 F1 = - 



32tt 2 



h 7b — 1 — ln47r 



- 2 J ( m )> 



and so the self energy in this scheme reads simply 

U F (p = 0,T) = to 2 



A TO 
In — 



32tt 2 



(19) 



(20) 



This is the zero temperature MS value, and we avoided the dangerous AT 2 terms again. In this case the condition of 
the constant physics reads: 



Am Ms 



ms 32tt 2 
which, at 0(A), implies 



h 7e — 1 — In 47r 

e 



h 7b — 1 — ln 4tt 



-^J(to f ) + 0(A 2 ), 



m MS = m F - 9 J ( m F)- 



This is another, also usual, form of the gap equation. The resummation again reads as 

'(p = 0,T)=n F (p = 0,T) 

This is a finite version of the thermal mass method. 



ytresummed { 
MS 



(21) 



(22) 



(23) 



C. Momentum dependent schemes 



What was said above is perfectly applicable to momentum independent resummations, as it was demonstrated in 
Ref. . To be able to treat momentum dependence we have to rethink what are really the conditions for a Lagrangian 
to be renormalizable. 

In fact, all theorems of renormalizability are based on the behavior of the Feynman diagrams when internal momenta 
are asymptotically large. One should prove that a diagram can be divergent if some of its internal momenta go to 
infinity on the same rate |l8j : we do not have new divergences when different momenta go to infinity at different 
rate. One proves that after removing the divergences of all sub-diagrams the diagram is superficially divergent with 
polynomial momentum dependence. Then the BPHZ scheme and forest formula go through, and we can perturbatively 
renormalize the theory [18}. 
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It is clear that all these theorems still hold, if propagators or vertices are momentum dependent, but in a way that 
the asymptotic behavior is the same as in the usual case. More precisely, the new momentum dependence should not 
cause any new divergences in any diagrams. 

For example, if we want to use momentum dependent "mass" term in <I> 4 theory, we have to take into account, that 
the most singular diagram contains one propagator and one integration (the tadpole). If the mass has asymptotic 
momentum dependence as 



m 2 (p) =m 2 R + 0{p-"<) , 
then the divergence of the tadpole diagram calculated with this "mass" reads 



d 4 p 



1 



d 4 p 



1 



(2tt) 4 p 2 - m 2 {p) J (2tt) 4 p 2 - m 2 R - O (p-~i) 



d 4 p 
(2^) 4 



O 



P 



(p 2 - m 2 R ) 2 



(24) 



(25) 



Here the second term yield finite result if 7 > 0; then it is not relevant for the issue of renormalizability. As a 
consequence we are allowed to use any momentum dependent mass term that approaches constant value polynomially 
at asymptotic large momenta. 

This means that we can define a consistent renormalization scheme, by splitting the bare mass term as 



ZrriQ — m 2 (p) + 8m 2 (p) 



(26) 



where m 2 (p) satisfies (I24|l . 

If the coefficient of the quadratic term is considered, ie. p 2 —m 2 , then the allowed modification is 0{p~ 2 ~ 1 ) relative 
to the leading p 2 term. A similar observation can be made for coefficients of higher powers. For example, if the 
quadratic coupling is denoted as A, then we can introduce a momentum dependent coupling as 



i Z 2 X bare (2 7 r) 4 J(Sp l ) $(pi)$(p 2 )$(p 3 )$(p4) = ^ (X(pi 



8X(j>i)) (2n) 4 8(X Pl ) S(pi)S(p2)$(p 3 )$(P4). 
If at high momenta the asymptotic behavior of \r is denoted as 

Hpi,P2,P3,P4) = A fl + 0((ma,xpi 
then the tadpole diagram is computed as 



d 4 p X(p,q,kJ) 



(2tt) 4 p 2 -m 2 {p) 



d 4 p 



A, 



O 



P" 



P 



-2-7 



P 



(27) 



(28) 



(29) 



where the second term is convergent if 7 > 0. 

So as a general conclusion we can state that a Lagrangian remains perturbatively renormalizable if the couplings 
are momentum dependent, provided they are suppressed by a factor 0{p~ 2 ~ 1 ) relative to the leading momentum 
power at high momenta. 

Finally we give two remarks. First, the order 0(p~ 2 ~ 1 ) claims nothing about logarithmic correction, since it has 
sub-power asymptotics. So, for example, an asymptotic behavior log(p)p~ 2 ~ 7 is still allowed. The second remark is 
that we may not stop at the level of the usual renormalizable operators. We can add ~ <J> 6 terms with coefficient of 
the leading term 0{p~ 2 ~ 1 ). This makes the renormalizable treatment of nPI diagrams possible. 



D. Momentum dependent renormalizable resummations 



Now we are ready to generally describe the resummations that are renormalizable at the same time. We restrict 
ourselves here to self-energy type resummations. 

We first define a scheme (M) that is particularly adequate for the problem at hand. In this case it is simply the 
choice of the momentum dependent mass counterterm Sm 2 I (p). Anticipating matching to MS scheme, we also use 
momentum dependent mass m 2 (p). To simplify the notation we introduce the renormalized free propagator in this 
scheme as 

G- A l (p)=p 2 -m 2 (p). (30) 

Now we perform perturbation theory in this scheme. Since it is a well defined scheme, all expectation values are finite. 
In particular the 1PI generating functional r^,f[<I>] is finite, too. This will depend on the propagator; to emphasize 
this dependence we shall write Fm[$; Gm.o]- 
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Now we translate the results to MS scheme. Since only the masses are modified (that is 8Z 2 and SX are the same 
as in MS) , the condition of the constant physics reads simply: 

TO Ms + 5m Ms = m A/0) + 5m 2 M {p). (31) 

In practical cases 5m 2 M (p) depends on the propagator Gm,o(p), and so this is a gap equation. Since the left hand side 
is momentum independent, m\[{p) should compensate the momentum dependence of 8m 2 M (p). If the counterterm 
has a polynomially vanishing momentum dependence, the mass will also have polynomially vanishing momentum 
dependence. In this case, according to the analysis of the previous subsection, the divergences in Sm 2 -^ and 8m 2 M (p) 
are the same. That means, that we have a finite relation between m 2 M (p) and rn 2 -^ - or, expressing in another way, 
the propagator in scheme M depends explicitly on the MS mass value: Gm,o(p', m Ms)' 

With this dependence all expectation values are the same in MS and in M schemes, if we calculated them at infinite 
order. At finite nth order they are the same only up to A™, the rest is the "resummation" as considered from the 
point of view of MS scheme: 



^ [<P\ = 1 m[9,Gm,q\ 

To make this, in essence very simple, strategy more tractable, we now work out some examples in the next sections. 



(32) 

Gm,o from constant physics 



III. RENORMALIZATION OF 2PI RESUMMATION 



A popular resummation method is the 2PI resummation; for a comprehensive description and further references 
cf. [a, • The physical idea is that the free theory should characterize the quasiparticle propagation in the best 
possible way, so we should use the exact propagator already in the free theory. Therefore in perturbation theory the 
internal lines are exact propagators, which also implies that it must not contain self-energy insertions. Since self- 
energy insertions are separated from the rest of the diagram by two lines, the above idea is based on diagrams that 
are two-particle irreducible. In the usual formalism it is achieved by introducing a generic propagator as an external 
quadratic current, and then perform a Legendre transformation on that. The requirement that the free propagator is 
exact can be formulated as a gap equation. 

There are different methods to renormalize 2PI resummation; these usually use BPHZ scheme, and rely on the 
renormalization of the zero temperature Bethe-Salpeter kernel 13]. The method proposed by this paper makes the 
2PI effective action finite in a different way, directly using the (generalized) counterterms of the original Lagrangian. 



A. The 2PI resummation as a renormalization scheme 



In the skeleton diagrams of the 2PI resummation the self-energy diagrams are missing. We can define a scheme, 
where this occurs, by appropriately choosing the mass counterterm. To understand the idea we recall the zero 
temperature on-mass shell scheme: here one cancels radiative mass corrections on the mass shell with help of the 
finite part of the mass counterterm, thus achieving that the tree level mass is the exact mass in the same time. The 
exact mass reads as 

m\ x =m 2 + ,5m 2 + n(fc)| fc2=mL , (33) 

and we have to ensure that it agrees with the free mass, ie. m 2 x = m 2 . That yields the following choice of the mass 
counterterm 

5m 2 + n(fc)| fe 2 =m 2 = 0, (34) 

therefore the renormalized self-energy correction is zero on the mass shell. To fully define the scheme, we should 
choose the finite parts of the other counterterms (<5A, SZ): we choose for them the MS values. 

We should generalize this scheme in a way that not only the mass, but the complete exact propagator is the same 
as the free propagator. We use a momentum dependent mass m(k), yielding the free propagator G _1 = k 2 — m 2 (k), 
together with momentum dependent mass counterterm 5m{k). From Schwinger-Dyson equation, using explicitly the 
mass counterterm, we can write 



G~J{k) = G~\k) - G- X (fc)(n(fc) + W(fc))G- 1 (fc). 



(35) 
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The condition that the free propagator is exact, G ex (k) — G(k), can be satisfied for all k if 

6m 2 (k) +n(fc) = 0. (36) 

This is in direct analogy with the on-mass-shell equation 13411 . We will call the scheme defined above as bare 2PI 
scheme. 

Two remarks are in order here. First we call the attention to the fact that II(fc) implicitly depends on the propagator 
G so on the momentum dependent mass m(fc), too. Secondly, if the bare 2PI scheme was a consistent scheme (as it 
is not, see later), the above condition should determine the finite part only, the regularized UV divergences should be 
the same as in any other schemes. 

If we want to have the result in MS scheme we must ensure that the 2PI and the MS schemes belong to the same 
class of constant physics. According to our previous analysis in Section III CI we should provide the same the bare 
mass in the two cases: 

Z-\m\k) +Sm 2 (k)) = ZdL(m^ + 6m^). (37) 

We can choose the same wave function renormalization counterterm in the two schemes. Then, using (J5SJ, we can 
write 

m 2 (k) - II[m, k] = + Sm^ (38) 
where we explicitly written out the mass dependence of the self-energy. Equivalently we can write 

m 2 (k)= m 2 ^ + [m, k] . (39) 

Here IIj^g[m, k] — II[to, k] — <^ TO ||gj a finite quantity (if the bare 2PI scheme was a consistent renormalization scheme). 
If we add a — k 2 to both sides we find 

G_1 ( fc ) =G M^- n Ms[ G ' fc ]- (4°) 

This equation is exactly the same as the 2PI gap equation that determines the self energy, or, equivalently, G(k). 



B. Improvement of the 2PI resummation 



Although the formalism of the 2PI resummation is appealing, it is inconsistent when using order by order. The 
self energy, namely, can yield corrections of type fc 2 lnfc, m 2 lnfc, which, if considered to be part of the propagator, 
leads to new, eventually temperature dependent, divergences in the calculation that cannot be canceled by usual 
counterterms. The related problem in the bare 2PI scheme is that the condition l|39|l requires that m 2 (k) contains 
momentum dependence ~ k 2 Ink or ~ m 2 Ink, which is not the allowed 0(fc~ 7 ) order. This means that the bare 2PI 
scheme is, in fact, not a consistent renormalization scheme. 

In the present framework we can accommodate to the requirement of renormalizability, if we employ 2PI resumma- 
tion in the IR regime, and leave the asymptotic momentum regime untouched. These schemes are already consistent, 
and we call them 2PI schemes. There are several possibilities to achieve this. 

The first is that we take into account the momentum dependence only below a scale M; above this scale we stay with 
the above discussed thermal mass approximation. That is we write for the mass counterterm, in high temperature 
approximation (cf. CH0 : 

f-n(fc) if|fc|<M 

Since at asymptotically high momenta we recover Lorentz invariance, we can use for |fc| the Euclidean norm. The 
condition of the constant physics (|39|l now reads as 



m 2 (k) = m2_+0(|fc| < M)IL m [m,k] + Q(\k\ > M) 



(42) 



This describes a scheme that provides 2PI resummation at low momenta, and super-daisy resummation at high 
momenta. Since we expect that the effects of the environment, and so the spectrum modifications are restricted to 
the IR regime, this choice hopefully catches the main physical points. 
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We can apply, instead of a sharp cutoff, a smeared cutoff function M(k) that is 1 for small momenta and vanishes 
asymptotically as fc~ 2 1 at least. Then we could write 

5m 2 (k) = -M(k)U(k) - (1 - M{k))- J(m 2 ), (43) 



and the condition of the constant physics reads 

m 2 (fc) = m^ + M{k)n m [m,k] + (1 - M(k)) 



— J(m 2 ) + 5m- 



MS 



(44) 



Another possibility, related to the technique applied by |l4j , is to subtract the problematic part in the asymptotically 
large momentum region: 

<5m 2 (fc) = -n[m, k] + U^ sympt [m R , k}. (45) 

The last symbol means the momentum dependence of the renormalized (in any scheme) vacuum self energy for 
asymptotically large momenta with mass ma; so this is an explicit function of uir and k. If condition Ij24(l is fulfilled, 
this equation ensures that the mass counterterm goes to a constant value at large momenta, since the momentum 
dependence in 5m 2 (k) must be slower than that of any functions present in Yl!^ mpt [mji, k]. The condition of the 
constant physics reads 

m 2 (fc) = mf^ + %s[m, k] np^V«, k\. (46) 

Here again l|24[l is a self-consistent Ansatz: if it is obeyed by m 2 {k), then the right hand side goes to a constant value 
at high momenta faster than any functions in HQ Svmpt [mji, k]. 

It must be emphasized that in any of the methods above we do not have complete 2PI resummation. To describe 
the high momentum regime we still have to rely on perturbation theory; but this regime is insensitive to matter effects 
and so any renormalization scheme applicable at zero temperatures is applicable here, too. 

We still have the task to characterize the asymptotically large momentum regime. Let us consider a two-point 
function in a renormalized theory and go to this regime. Here the only dimensionful quantity is the momentum, 
which is made dimensionless with help of the renormalization scale. So we can write the self energy as 

S (asym)( fc ) _ k 2 Z ^^{k / [l) + W 2 Z^ Sym \k / . (47) 

Therefore the dangerous terms in the asymptotic regime can be obtained from the scale dependence of the self- 
energy. This means that we should count with logarithmic dependence, and the coefficients of the logarithms are the 
renormalization group beta-functions. 

IV. 2PI SCHEME IN SIX DIMENSIONAL $ 3 THEORY 

The simplest theory where we can demonstrate how the 2PI scheme introduced above works, is the 6D $ 3 theory 
at finite temperature. The Lagrangian of the theory is written as 

£ = -(d$) 2 - — $ 2 - -$ 3 + counterterms. (48) 

Although the theory is not well defined in the path integral sense, it behaves well in the perturbation theory, and, 
moreover, it shows considerable resemblance with the perturbation theory of the gauge models. We will use one loop 
level resummation; we will have momentum dependence already at this level, so a momentum dependent resummation 
is needed. On the other hand the divergence structure is rather simple, which makes the treatment easy. 

The most important peculiarity of the present treatment is that we allow momentum dependence in the mass 
counterterm, which implies momentum dependent mass. So the free Lagrangian can be written in momentum space 

C = i(fc 2 - m 2 (fc))$(fc)$(-fc) - -5m 2 (k)$(k)$(-k), (49) 



in such a way that the bare mass is untouched 

m bare = 7/1 \^ ) T V" ) - '""MS U/,t MS 



m 2 (fc) + 5m 2 (k) = + 8mL^. (50) 
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We will use real time formalism in R/A basis W]- Here the original fields ($1, $2) are replaced by (<i> r ,$ a ) via 
the definition 

$1 = *r + ^a, $2 = $r - ^a- (51) 

The propagator is a 2 x 2 matrix, where, in this basis G aa = for the exact Green's function, and also G* a (k) = G ar (k) 
is true in the Fourier space. The interaction Lagrangian reads 

-£,= ! [<&?-#!] =|$2$ a + A$3 (52) 

In the momentum dependent resummation we should use the most generic possible mass matrix. Since we are in 
equilibrium, all the propagators can be derived from the retarded propagator. This implies relations between the 
generic self energies, which we should respect when we choose the mass matrix. As it was pointed out in 20], and 
summarized in Appendix lAl the sclf-cncrgics satisfy 

Re En = — Rc E22 = Re E.R, Im En = Im E22 = (1 + 2n R ) Im 

£ 12 = — 2msImEfl, E 2 i = — 2i(l + n B ) ImEjj. (53) 

Here we used the notation E^ = E ar which is the retarded self energy. This is because it parameterizes the retarded 
Green's function as 

G- R 1 {k) = G- R ] {k)-^ R {k). (54) 

Choosing the mass matrix according to l|53[) . the only independent parameter is the m R (k) retarded mass, and the 
free retarded propagator reads G RQ (k) = k 2 — m R (k). This implies for the free spectral function 

—2 Im m 2 f k) 

go (k) = DiscGnAk) = -2hnG R , {k) = (fc2 _ Rem2#))2 | (Imro 2 (fc))a ■ ( 55 ) 

From the properties of Qo(k) coming from causality, this equation requires that lmm R (k) must be an odd function of 
ko, and also lmm R (k) < must be true for ko > 0. In the MS scheme it is satisfied by choosing lmm R (k) = — e sgnfco- 
Since Gr^ should be causal, m R must be either local or causal, ie. Kramers-Kronig relation holds true for it: 

oc 00 

m»(fa,k)= / ^' Imm ^' k) = /^M-^^y, (56 ) 
J 7r k - lj + is J it (k + ie) 2 - uj 2 

-00 

Therefore Im m 2 R and the local parts of m 2 R will completely determine Gr^ . 
For the counterterms we know that the infinite part is fixed, namely 

C cttdiv = \{k 2 5Z 2 8ml v ){*{ - <S> 2 ) - - d> 2 ) 2 . (57) 

This implies KeSm R = Sm 2 iv + finite, and \m.8m 2 R = finite. We choose for 8Z 2 and Sg their MS values. 

If we do perturbation theory in a generic scheme, all the diagrams are expressed through the propagators. Since 
all the propagators are functionals of g, thus every diagram is a functional of g. In particular, the one loop retarded 
self-energy in <I> 3 model can be written as 

Efl(fc) = g 2 I R (k) + Sm 2 R (k) + k 2 SZ 2 , (58) 

where / is the "bubble" diagram 

U R (k) = J j—g iGrrfiip) iGra : o( k - (59) 

At finite temperature the KMS relation implies 

iG rrfi {p) = Q + n(po)^j g (p). (60) 
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Thus the imaginary part of the retarded self-energy reads 



lmI R (k) = ~J So(p)Qo{k-p), 



(61) 



indeed a functional of g alone. 

For later convenience, this contribution can be divided into three parts, a "regular" , a "singular" and a "divergent^ 
part: 



i R (k) = it(k) + 4 DS (k) + i^(k). 



(62) 



The divergent piece is defined that in a certain regularization this part goes to infinity as the regularization parameter 
vanishes. The singular part is defined as a contribution that grows at large momenta. Since the bubble contains 
no divergent sub-diagrams, these contributions can be calculated from the zero temperature result. In dimensional 
regularization this reads as 



t t=o 
l r 



(k) 



1 



2(4tt) 



- t J t In 4tt j ( -y + in 



A I 1.2 ~C\ \ -L 2 \\ ~k 2 x{l -x) + m 2 
ax {—k x(l — x) + m ) In 

A* 



fco — *ko-\-is 

(63) 

In the real part we have to take the modulus of the argument of the logarithm. The integral can be evaluated, the 
result reads 



where 



ReI R ^(k) = 



1 



2(4tt) 



1 8 , m 2 
- - IE + - - In j 



k 2 



k 2 T 3 



h m z } 1 In 



i + r 



i -r 



Im ir°W = -6(fc 2 - Am 2 ) 

R y ' y ' 128tt 2 6 ' 



(64) 



r = \ l- 



4m 2 



(65) 



What is particularly interesting is its large fco behavior; for a fixed k it can be obtained as the large k behavior: 



Re/Jr°(fc -> oo) 



1 



2(4tt) s 



1 8 , k 2 

- - IE + - - In 2 

e 3 Att^ 



k 2 



ImJ^-°(fc -» oo) 



sgn(fco) / fc 



(i 



1287T 2 

For the divergent pieces we choose the following expression 

g 2 I R iv (k) = -k 2 8Z 2 - m 2 8Z- 

where 



&7 2 - 1 
6Za "I2(47r) 3 

As a singular part we choose: 



1 8 , „ 

- - IE + ~ + hi 47T 



<5^1 = - 



1 



2(4tt) 3 



1 13 , „ 
7b + — + In 47r 

e 6 



(66) 



(67) 



(68) 



k 2 



2(4tt) 3 V 6 



In 



M 2 



fco — >ko-\-ie 



(69) 



where M 2 is an arbitrary scale. 

The regular piece, from (|62|) . can be obtained as = I R — I R 1V — I R ne - Its imaginary part, from the definition 
(|59|l and from the choice of the singular part of the diagram, reads as 

1 f d 6 p (\ 



■x + n iPa) Qo(p) Qo(k -p) + 



-r^Q(k 2 - M 2 ) f—-m 
2(4tt) 3 V ' \ 6 



(70) 
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From (|59|l and 161(1 one can see that the real part can be obtained from the Kramers-Kronig relation 

— oo 

This integral is convergent, since lmI R s (k) ~ 1/fc 2 for large fco values. 
The self-energy reads with these choices as 

Z R (k) = -k 2 5Z 2 - m 2 SZ 2 l + g 2 I R hlg (k) + g 2 I R & (k) + Sm 2 R (k) + k 2 SZ 2 . (72) 

In the 2PI scheme we choose the countcrterms in a specific way, namely we choose 

Sm 2 R (k) = m 2 SZ 2 m - g 2 I R s (k), SZ 2 = SZ 2 . (73) 

Thus almost all terms are canceled in the self energy; what remains is 

where the Landau prescription (Aft — > fco + is) is implicitly understood in the formula. We can observe that this 
expression is almost constant in the IR regime where k 2 <c M 2 . For convenience, in the numerical calculations we 
have chosen /i = M, and so this constant is ~ ln(l — k 2 /M 2 ) <C 1 for k 2 -C M 2 . 
The full retarded propagator in this scheme reads as 

n 2 (k 2 \ M 2 - k 2 
G R \k) = k 2 - ml(k) - ^ [L - m j In ^_JL, ( 75 ) 

Our reference scheme will be the renormalization scheme defined by the mass counterterm Sm 2 r = m 2 6Z 2 and 

rGl ft i 

wave function counterterm SZ 2 c{ = SZ 2 ,. For comparability of the results we have to ensure that the bare Lagrangian 
is the same in the resummed and in the reference scheme. Since the wave function renormalization is the same in the 
two schemes, we have to ensure that 

Z 2 m 2 aic = m 2 ci + 6m 2 ci = m 2 R (k) + Sm 2 R (k) = m 2 R (k) + m 2 5Z 2 m - g 2 q s (k). (76) 

This equation implies 

m R (k) = m 2 Ief + g 2 I^(k). (77) 

This is a finite gap equation that can be solved numerically. 

To determine the value of m 2 cl we used the renormalization condition that the exact mass should be a predefined 
value m 2 at zero temperature; ie. we required 

m R (k 2 = m 2 ,T = 0) = to 2 . (78) 

From l|77Jl we can read off the value of m 2 o{ , which later can be used at finite temperature calculations. 

Before giving the results of this calculation, we mention that formally it is also possible to choose /i 2 = M 2 — k 2 . 
With this choice the 2PI scheme is exact. The validity of the above formulae can be maintained, if we choose the 
same scale in the reference scheme. Then the condition of the constant physics is modified as: 

m 2 R (k)=m 2 M+g 2 (n)r^(k). (79) 

Here TO 2 cf (/i) and g 2 (n) can be explicitly calculated using the standard renormalization group. This choice may be 
advantageous in case of asymptotically free theories. In the $ 3 model, however, we have to be careful because of the 
UV Landau pole. 



A. Results 



For numerical calculation we used the spatial rotational invariance of the propagators, and so at finite temperature 
they are two-dimensional functions, with variables po and |p|. The integrations are performed on 2D lattices, which 
are mapped with a continuous function to the infinite 2D space. 
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The main result of the calculation is the retarded self-energy. Its imaginary and real part at zero and at finite 
temperature can be seen on Fig. ^ I n this run the parameters are g = 20, M/m — 10 and for the finite temperature 
part T/m = 1. The first plot shows the imaginary part. We can see that at zero temperature (curve (a)) it is zero 
below the 2-particle threshold ko/m = 2. With dotted line we plotted the one loop result from Q64[l. We can see that 
below the 3-particle threshold (ko/m = 3) the agreement is rather good. The most prominent phenomena of this plot 
is that at finite temperature the spectral function is non-vanishing everywhere (except zero) - this is what we expect 
on general arguments. One can observe the smeared-out remnant of the Landau-damping (for ko < k), but between 
the Landau-damping region and the zero temperature threshold the spectral function is still large. Remarkably, in 
this regime the imaginary part is almost constant. 

On the second plot of Fig. 2] the real part of the full self energy from the same calculation is depicted. We can see 
that at zero temperature it is correctly renormalized, we indeed obtain T,(ko/m = 1) = to 2 . At finite temperature 
we obtained larger real part: this means that the temperature dependent mass correction is positive. In the thermal 
mass approximation the finite temperature curve would be above the zero temperature part by a constant shift. But, 
according to this figure, there is considerable difference between the fco — contribution (Debye-mass) and the ko = to 
contribution (on-shell mass correction). 

On Fig. |2 the temperature dependence of the quasiparticle properties can be seen. In the first plot the thermal 
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FIG. 2: The first plot describes the thermal mass correction. The dashed curve (a) corresponds to the Debye mass (ie. ko = 
case), the solid curve (b) is the on-shell value. In the second plot the on-shell imaginary part of the self energy can be seen. 



mass correction is shown. The two curves correspond to the thermal mass corrections evaluated at fco = (curve 
(a): this is the Debye- or screening mass) and on mass shell (curve (b): this is the quasiparticle mass). There is a 
rather big difference between the two values, which indicates that a single effective mass term cannot be a satisfactory 
description. In the regime T/m > 0(1) both corrections can be well fitted by a quadratic function of the temperature 
(aT 2 + 6), but in the small temperature regime the curve is much shallower. This means that although the formula 
TO 2 ff ~ to 2 + aT 2 is valid at high temperatures, but to is smaller than the zero temperature mass. If this is a robust 
phenomenon persisting to models with spontaneous symmetry breaking, it should lead to increasing of the critical 
temperature of a phase transition as compared to the pure perturbation theory. In the current case (g = 20) the mass 
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modification is about 5%. 

In the second plot of Fig. [21 the on-shcll value of the imaginary part of the self energy can be seen as a function of 
the temperature. This is purely resummation effect: at one loop level there is no on-shell damping. The result again 
fits well to a quadratic temperature curve (aT 2 + b) in the regime T/m > 0(1), and it is shallower than quadratic 
for small temperatures. At high temperatures it leads to a quasiparticle damping which depends linearly on the 
temperature. 



When we leave the realm of the strict perturbative calculations and try to implement exact, or approximate 
resummation formulae (like Schwinger-Dyson equations or 2PI resummation equations), we encounter new type of 
divergences in the calculation. We have then two ways to proceed. One is that we try to find out the structure 
of the new type of divergences, then resum the counterterms in the same way, thus rendering the calculation finite. 
An alternative way is to realize that divergences come solely from the ultraviolet regime which is perturbative (if 
it is perturbative at zero temperature), while the interesting results are at low energies. This allows to change the 
resummation equations in a way that they do nothing at high energies, but do resummation in the physically important 
IR regime. 

In this alternative way the divergence structure does not vary with the environment, which means that, in fact, 
we use a special renormalization scheme. Technically this scheme works in a way that we choose the finite parts 
of the counterterms (which are free to choose) to cancel the IR part of the IR sensitive diagrams. The remaining 
perturbation theory is no more IR sensitive, so we can safely use it. The results are, however, obtained in a weird 
scheme which depends on the environment, while we would need the result in a well defined, fixed scheme (eg. in 
MS). But results calculated in two schemes can always be connected by renormalization group transformation: we 
only have to ensure that the original bare Lagrangian is the same in the two cases. To satisfy this requirement we 
have to write up constraint equations which are, in the original resummation language, the gap equations. 

This strategy can be used for momentum dependent schemes, too. In particular one can derive a 2PI scheme where, 
in the IR regime, all the two particle reducible diagrams are canceled in the perturbation theory. The requirement 
of keeping the bare Lagrangian constant can be translated to a consistency equation for the momentum dependent 
mass. If we treat the UV characteristics of the theory well, these equations provide finite solutions. 

To demonstrate in an example how this strategy works, we computed the 2PI resummed renormalized self energy in 
the six dimensional $ 3 model. We used MS scheme as reference scheme, and we studied finite temperature dependence 
of different physical quantities. The model at relatively large coupling (g = 20) still can be characterized by finite 
lifetime quasiparticles, but the spectral function (density of states) is nonzero for all nonzero frequencies. Interesting 
result is that the on-shell thermal mass and the zero frequency (Debye) mass are quite different, in the thermal mass 
correction there is almost a factor of 2 between the two values. 

To pursue this project there are two ways open. Since we worked out the technicalities, we will be able to use 
renormalized 2PI analysis for other models, too. On the other hand, with the same logic we can discuss the 4PI or 
even higher point irreducible approximations in a renormalized way. 



The author is grateful to A. Arrizabalaga, T.S. Biro, A. Patkos, J. Polonyi, U. Reinosa and Zs. Szep for useful 
discussions. This work was supported by Hungarian grant OTKA F043465. 



Most easily we can arrive to the relations between the self-energies in the following way. First we determine the 
relation between the self energies in the CTP and R/A formalism. For this we write up the quadratic part of the 
effective Lagrangian in equilibrium: 



V. 



CONCLUSIONS 
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APPENDIX A: THE MASS MATRIX AT FINITE TEMPERATURE 
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To change to the R/A formalism we use 

$ (i,2) = q>r ± i $ a Thcn the effective La 

grangian can be written as 

A2) _ i($ r , $ a ) ( k 2 - E ra \ (9 

where 



^'-2 ' V-E. E- *• ■ (A2) 



En 


+ Ei 2 


+ E 2 i 


+ E22 


— Z-if<y — 
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+ E 12 


— E21 
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— E12 


+ E21 


— E22 


— 2S ra 


En 


— E12 


— E21 


+ E22 


= 4S Qa , 



(A3) 



where we used the information that G aa = 0, and 


so S^.^ - 


= 0. The ab< 


E fl 


= E Qr = 


= E U - 


f S12 


Ea 


— E ra - 


= En - 


f E 2 i 


Ek 


= E aa = 


E21 


+ E12 
2 



(A4) 



On the other hand, inverting the kernel in the R/A formalism yields the propagators 



11S 

Gr = T2 v^> Ga = U2 {Gk = v — ( A5 ) 

- E_r k 2 - E A S K - S A 

where p = iGn — iGa- Since GJj = Ga, therefore S|j = E^. Moreover, since zG/f = (i + n)g, so S oa = (| + 
n)(— 2i IrnSj}). That means 

En = ^ SA - E aa = Re E fl + i(l + 2n) Im E fi 
E 22 = - I + E aa = - Rc S fl + i(l + 2n) Im S fl 
E12 = h E aa = -2mImE iI 

E 2 i = - SK ~ SA + S aa = -2»(1 + n) Im S fl . (A6) 
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